Fermionic mode entanglement in quantum information 
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We analyze lermionic modes as fundamental entities for quantum information processing. To this 
end we construct a density operator formalism on the underlying Fock space and demonstrate how 
it can be naturally and unambiguously equipped with a notion of subsystems in the absence of a 
global tensor product structure. We argue that any apparent similarities between fermionic modes 
and qubits are superficial and can only be applied in limited situations. In particular, we discuss 
the ambiguities that arise from different treatments of this subject. Our results are independent of 
the specific context of the fermionic fields as long as the canonical anti-commutation relations are 
satisfied, e.g., in relativistic quantum fields, or fermionic trapped ions. 



I. INTRODUCTION 

Fermionic systems have been analyzed as agents for 
quantum information processing in a multitude of stud- 
ies, ranging from discussions of fermionic modes of 
relativistic quantum fields [1-11], over fermionic lat- 
tices [12], and fermionic Gaussian states [13], to dis- 
cussions of the entanglement between fixed numbers of 
indistinguishable particles [14-24]. In the latter case, 
only pure states of fixed particle numbers are consid- 
ered and a selection of entanglement measures are avail- 
able, see, e.g., Ref. [20]. However, these restrictions 
seem to be much more limiting than required. From 
the point of view of quantum information theory it is 
natural to ask for an extension to incoherent mixtures 
of quantum states. Furthermore, from the perspective 
of a relativistic description particle numbers are not 
usually conserved, i.e., the particle content of a given 
pure state is observer dependent [25]. Only the coher- 
ent superpositions of pure states of different (electric) 
charges are forbidden [26]. 

In the light of this fact it is therefore reasonable to con- 
sider the entanglement between fermionic modes, in a 
similar way as this is conventionally done for bosonic 
modes, e.g., for Gaussian states [27]. However, due to 
the Pauli exclusion principle fermionic modes are nat- 
urally restricted to two degrees of freedom, i.e., each 
mode can be unoccupied or contain a single excita- 
tion. This has provided many researchers with an ad 
hoc justification for the comparison with qubits - two 
level systems used in quantum information. The central 
question then simply is: Can fermionic modes be con- 
sidered as qubits? The short answer to this question is 
"No" , but in limited situations certain techniques from 
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the study of qubits can be applied in fermionic sys- 
tems. This has incited debates among scientists, see, 
e.g., the exchange in Refs. [28-30]. In the following we 
shall give a more precise answer to the question above, 
along with a detailed description of the problem. The 
supcrselection rules further restrict the possible opera- 
tions that can be performed on single-mode subsystems 
and it was argued that this should lead to a modified 
definition of the entanglement between modes [20]. At 
least for fixed particle content this problem can be cir- 
cumvent [31]. Moreover, even if quantum correlations 
are not directly accessible, a transfer of the entangle- 
ment to systems that are not encumbered by such re- 
strictions should be possible, thus justifying the use of 
unmodified measures for mode entanglement. 
The main aim of this paper is to establish a clear frame- 
work for the implementation of fermionic field modes as 
vessels for quantum information tasks. To this end we 
present an analysis of the problem at hand, i.e., how 
the modes in a fermionic Fock space can be utilized as 
subsystems for quantum information processing. We 
present a framework that is based on simple physi- 
cal requirements in which this can be achieved. We 
further discuss the issues and restrictions in mapping 
fermionic modes to qubits and we show how previous 
work and proposed solutions, e.g., invoking superselec- 
tion rules [29], fit into this framework. 
The article is structured as follows: We start with a 
brief discussion of the description of fermionic Fock 
spaces in Sec. II and how density operators arc con- 
structed on such spaces in Sec. III. We then go on 
to formulate the "fermionic ambiguity" that has been 
pointed out in Ref. [28] in Sec. IV. Subsequently, we 
reinterpret this as an ambiguity in the definition of 
mode-subsystems, which can be resolved by physical 
consistency conditions, in Sec. V. Finally, we discuss 
the implications for the quantification of entanglement 
between two fermionic modes in Sec. VI, before we in- 
vestigate situations beyond two modes in Sec. VII. 
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II. THE FERMIONIC FOCK SPACE 

Let us consider a (discrete) set of solutions ip n to a 
(relativistic) field equation, e.g., the Dirac equation. In 
the decomposition of the quantum field ip each mode 
function ip n is assigned an annihilation operator b n and 
a creation operator b^ , which satisfy the canonical anti- 
commutation relations 

{b m ,b\ l }= 6 mn , (la) 

{b m ,b n }= {bl,bU = 0, (lb) 

where { . , . } denotes the anti-commutator. One may 
introduce a different notation for the annihilation and 
creation operators for modes with opposite charge, see, 
e.g., Ref. [8], but for the purpose of our analysis here 
this is inconsequential and we can work only with b n 
and b n . The creation operators b^ n , acting upon the 
vacuum state | ), will populate the vacuum with a 
single excitation, i.e., 

\i>n)=bl\0} , (2) 

while the vacuum is annihilated by all b n , i.e., b n | ) = 
Vn. As can be quickly seen from this property and 
Eq. (la), the states \ip n ) are orthonormal. The states 
| ip n ) further form a complete basis of the single particle 
Hilbert space Hi- P , whereas ) E Ho- p ^ Hi- p . A 
general state in H.\- p has the form 

\^ P ) = X> , (3) 

i 

with £\ \m | 2 = 1 such that ( V 1_p | ip X ~ p > = 1. Let 
us now turn to states of multiple fermions. A second 
fcrmion can be added to the state (2) by the action of 
another creation operator b m , i.e., 

blfk |0)<X \lP m M ■ (4) 

Clearly, the anticommutation relations (1) require the 
two-fermion state to be antisymmetric with respect to 
the exchange of the mode labels m and n. We therefore 
define 

\lpm,1pn) = blbl |0) = |Vm)A|VO (5) 
= ^ (| 1pm ) ® | Ipn) ~ Ipn ) ® | 1pm )J . 

The two-fermion states are thus elements of the anti- 
symmetrized tensor product space of two single-fermion 
Hilbert spaces, i.e., 

H 2 -p= s(Hi- p ®Hi- p ) t (6) 

and a general state within this space can be written as 

|V 2 " P )= X>j , (7) 

i,3 



where the coefficients /iy form an antisymmetric ma- 
trix. States with more than two fermions can then 
be constructed by anti-symmetrizing over the corre- 
sponding number of singlc-fcrmion states. Finally, the 
n-mode fermionic Fock space T n is simply given as 
the direct sum over all fermion numbers of the anti- 
symmetrized Hilbert spaces, i.e., 

n 

Fn{U x . p ) = S(UT P ) (8) 

rn—l 

= Hq- p © Ui~ p © S(Ui-p ® Hi-p) (S ... , 

where H® m denotes the m-fold tensor product and we 
write Ho- P as Hf_ p . A general state in the space T n 
can be written as 

n 

1***)= A*o[ 0> ©£w|^i> (9) 

i=l 

8 72^2^k\ipj) A\*p k ) © ... . 

Let us now simplify the notation. From now on we will 
denote states in the fermionic Fock space by double- 
lined Dirac notation, i.e., || . }) instead of | . }, where 
the antisymmetric "wedge" product is implied when 
two vectors are multiplied, i.e, || ■ }} || ■ }} = || • )} A || . }). 
Furthermore, let us use the common "occupation num- 
ber" notation and write 1„ instead of ip n to denote an 
excitation in the mode n. Finally, we omit the sym- 
bol for the direct sum and simply keep in mind that 
states with different numbers of excitations occupy dif- 
ferent sectors of the fermionic Fock space. With this 
convention in mind we can rewrite Eq. (9) as 

l*»=Mo|0)}+X>|l*}) (10) 

i=l 

+ 7il>;*|l;>>|l*>> + --- • 

For the adjoint space we use the convention (compare 
to Eq. (5)) 

«1„|<(U| := (0\\b n b m = (blb n I0))) 1 (11) 

= - ^ ( ( Ipn | ® ( 1pm I - ( 1pm I ® ( Ipn |) , 
which allows us to write 

((lm|| (( In I li )) ||lj» = 8 m S mj - S nj S mi . (12) 

This notation is more convenient for computations in 
the fermionic Fock space. It should be noted that, in 
standard quantum information notation, the position 
of a "ket" corresponds to a particular ordering of the 
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subspaces with respect to the tensor product structure 
of the total space. Here, however, there is no tensor 
product structure corresponding to different modes ac- 
cording to which the vectors || . )) can be naturally or- 
dered. 



III. DENSITY OPERATORS IN THE 
FERMIONIC FOCK SPACE 

In complete analogy to the usual case of mixed states 
on tensor product spaces we can now construct incoher- 
ent mixtures of pure state in a fermionic Fock space. 
Let us first consider the projector on the state || "J }} 
from Eq. (10), i.e., 

I*M*I= ImoI 2 II o»«o || + £>,4 \uKij\ 

j,j',k,k' 

+ E(^s ii 1^ » « o ii + h.c.) + ... . 

(13) 

Let us check that such an object satisfies the criteria 
for a density operator: 

(i) It can be immediately noticed that (13) provides 
a hermitean operator. 

(ii) The normalization, i.e., Tr(| ^ ))(( * ||) = 1, is 
guaranteed by the normalization of || \& )) . In other 
words, the trace of (13) is well defined and in- 
dependent of the chosen (complete, orthonormal) 
basis in J- . 

(iii) Positivity: Finally, the eigenvalues of || \1/ ))(( \I> || 
are well defined, i.e., (13) can be represented as a 
diagonal matrix with diagonal entries {1,0,0,...}, 
which clearly is a positive semi-definite spectrum. 

We can then simply form incoherent mixtures of such 
pure states using convex sums, i.e., 

e=X>»l*»M*»l (14) 

n 

where X)nP" = 1, to construct the elements of the 
Hilbert- Schmidt space J-Ls{J~) over the fermionic Fock 
space. Properties (i) and (ii) can trivially be seen to be 
satisfied for such mixed states. The positivity of (14), 
however, requires some additional comments. The op- 
erator g can be diagonalized by a unitary transforma- 
tion U on which in turn can be constructed from ex- 
ponentiation of hermitean or anti-hcrmitcan operators 
formed from algebra elements b n and . Operationally 
this procedure is rather elaborate. A simpler approach 
is the diagonalization a matrix representation of g. As 



we shall see in Sec. IV, the matrix representation of 
g is not unique, but all possible representations tt^q) 
are unitarily equivalent, such that their eigenvalues all 
coincide with those of g, i.e., 

spectr (71^(0)) = spectr(g) Mi. (15) 

IV. THE FERMIONIC AMBIGUITY 

Let us now turn to the apparent ambiguity in such 
fermionic systems when quantum information tasks are 
considered. It was pointed out in Rcf. [28] that the 
anti-commutation relations (1) do not suggest a natural 
choice for the basis vectors of the fermionic Fock space 
for the multi-particle sector, i.e., for two fermions in 
the modes m and n, either 

|l m »|l„» Or |l4||l m }} = -|lr4l|ln>> (16) 

could be used to represent the physical state. This 
becomes of importance when we try to map the states 
in a fermionic n-mode Fock space to vectors in a ?i-fold 
tensor product space, i.e., 

TTi : T n — > Ui ® . . . <g> U n (17a) 
W) ^ I ) (17b) 
g ^ TTi(g) (17c) 

where the spaces Hi = C 2 (i = l,...,n) are identi- 
cal, singlc-qubit Hilbert spaces. The mappings TTi are 
unitary, i.e., (( <j>\ ip }) = ( <j> (l) | ) an d Tt(ga) = 
Tr(7ri(p)7ri(cr)). This implies that the maps 7Tj for dif- 
ferent i are unitarily equivalent. In particular, the dif- 
ferent matrix representations TTi(g) are related by mul- 
tiplication of selected rows and columns of the matrix 
by (-1). 

In the language of quantum information theory the 
states ip^i) arc related by global unitary transforma- 
tions. It thus becomes apparent that the entanglement 
of TTi(g) with respect to a bipartition 

n llm \u ptm+1 g)...® (is) 

will generally depend on the chosen mapping. 
Clearly, this is an unfavorable situation, but the in- 
equivalence of entanglement measures for different such 
mappings has been noted before, sec, e.g., Rcfs. [13, 23, 
29], while other investigations [8-10] did not suffer from 
any problems due to this ambiguity. Recently, the au- 
thors of Ref. [29] suggested that the ambiguity can be 
resolved by restrictions imposed by charge superselec- 
tion rules, while Ref. [11] suggested a solution by en- 
forcing a particular operator ordering. We will discuss 
both of these approaches in Sec. V, where we present 
simple and physically intuitive criteria for quantum in- 
formation processing on a fermionic Fock space. Most 
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importantly, we will show in Sees. V and VII that map- 
pings of the type of (17) can only be considered to be 
consistent when limiting the analysis to two fermionic 
modes obeying charge supersclcction, but not beyond 
this regime. 



V. THE PARTIAL TRACE AMBIGUITY 

While the sign ambiguity in the sense of the different 
mappings 7Tj is the superficial cause of the issue we 
want to discuss now a separate, and in some sense 
more fundamental problem: partial traces over "mode 
subspaces" . We arc interested in the entanglement 
between modes of a fermionic quantum field. How- 
ever, in the structure of the Fock space, there is no 
tensor product decomposition into Hilbcrt spaces for 
particular modes, see, e.g., Eq. (5). Only a tensor 
product structure with respect to individual fermions is 
available, but, since the particles are indistinguishable 
the entanglement between two particles in this sense 
has to be defined very carefully [20]. This issue is not 
unique for fermions and is sometimes rcfereed to as 
"fluffy bunny" entanglement, see Ref. [21]. 

For the decomposition into different modes we only 
have a wedge product structure available. In Ref. [20] 
the authors suggest that entanglement should be con- 
sidered with respect to this special case of the "braided 
tensor product" . As far as the construction of the den- 
sity operators with respect to such a structure is con- 
cerned we agree with this view, see Sec. Ill, and no 
ambiguities arise regarding the description of the total 
n-mode system. However, the crucial problem lies in 
the definition of the partial tracing over a subset of the 
n modes. This is best illustrated for a simple example: 
Consider a system of two fermionic modes labelled n 
and k'. A general, mixed state of these two modes can 
be written as 

Q K * = «i II » « || + a 2 |ltfMM (19) 

+ a 3 ||U«U + <* 4 \\u))\\U>KU'\\{U\\ 
+ (A I OHM + h | OHM 
+ ft l|o))«M|«i K | + fh ju.HU 
+ (3 5 wwHU'WiUW + & IIU«i,'ll«U 

+ h.c.\ 

where appropriate restrictions on the coefficients on € 
K and f3j £ C apply, to ensure the positivity and nor- 
malization of q kk i. Here we have, for now, disregarded 
supersclcction rules. Let us now determine the corre- 
sponding reduced density operators (on the Fock space) 
for the individual modes k and k' . Usually one would 



select a basis of the subsystem that is being traced 
over, e.g., for tracing over mode k' one could choose 
{|| )) , 1 1 K ' ))}. This clearly cannot work since basis vec- 
tors with different number of excitations are orthogo- 
nal. We thus have to define the partial trace in a differ- 
ent way. This is equally true for bosonic fields as well. 
However, in contrast to the fermionic case, no ambigu- 
ities arise in such a redefinition for bosonic fields. For 
the diagonal elements of the reduced fermionic states 
the redefinition of the partial trace is straightforward 
as well. These elements are obtained from 

Tr m (||0))«0|):= |0»«0||, (20a) 

Tr ro (I l n )) (( l n |) := 8 mn I U » (( U I (20b) 

+ (1-<W>||0}}<(0|| , 

Tr m (||l m ))||l„))((l„||((l m ||) := |l„Mln| (jn ? n) , 

(20c) 

where n,m = k,k'. While the diagonal elements are 
unproblematic and do not suffer from any ambiguities, 
we have to be more careful with the off-diagonal ele- 
ments. Three of these will not contribute, i.e., 

Tr m (|l m ))((l„|) = Tr TO (|| 0»((1 TO || ((1„||) (21) 
= Tr m (|l„))((l m |((l„|) = 0, 
and two more are unproblematic as well, i.e., 

Tr m (|| »((l n |) := (1 - S mn ) || »«1„ || . (22) 
The last element, 

Tr m (||l m )}((l m ||((l n ||) = -Tr m (||l m ))((l n ||((l m ||) 

= ± | Mini , (23) 

however, presents an ambiguity. If a mapping 7Tj to 
a two-qubit Hilbert space is performed, the choice of 
map will determine the corresponding sign in the par- 
tial trace over either of the qubits. The differences in 
entanglement related to the fact that TVi(g) and irj(g) 
are related by a global unitary are thus explained by 
the relative sign between the contributions of Eq. (22) 
and Eq. (23) to the same element of the reduced den- 
sity matrix. 

However, simple physical requirements restrict the 
choice in this relative sign. Any reduced state formal- 
ism has to satisfy the simple criterion that the reduced 
density operator contains all the information about the 
subsystem that can be obtained from the global state 
when measurements are performed only on the respec- 
tive subsystem alone. 

Let us put this statement in more mathematical terms. 
For any bipartition A\B of a Hilbcrt space H (with re- 
spect to any braided tensor product structure on H) 
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and any state p E % the partial trace operation Tr^ 
must satisfy 

(O n (A)) p = (O n (A)) Tl . B{p) , (24) 

where (0) p denotes the expectation value of the op- 
erator O in the state p and {O n (A)} is the set of all 
(hermitean) operators that act on the subspace A only. 
For the operator g KKi i from Eq. (19) the condition (24) 
can be written as 

Tr(O n ( K )g KK >) = Tr(O n (K)g K ) , (25) 

where g K = Tiv(g KK <). This consistency condition 
uniquely determines the relative signs between different 
contributions to the same elements of g K . Let us con- 
sider the (hermitean) operators (b K + bl) and i{b K — b\). 
Their expectation values for the global state g KK > arc 
given by 

Tr( (b K + b\)g KK , ) = 2 Re( (3 2 + ft ) , (26a) 
Tr( i(b K - b\)g Kti , ) = 2 Im( ft + ft ) . (26b) 
For the mode k! on the other hand we compute 

Tr( (b K > + b[,)g KK , ) = 2Re(ft - ft ) , (27a) 
Tr(i(b K , - b[,)g KK , ) = 2Im(ft - ft) . (27b) 

Equations (26) and (27) determine the sign in Eq. (23) 
and we find the reduced states 

g K = Ti V (q kk >) = (oi + on) II ))(( | (28a) 
+ (a 3 + « 4 ) |1«M1«| 
+ ((ft + ft) || » «1« 1| + h.c.) 

g K i = Tr K (g KK ,) = {a x + a 3 ) || ))(( || (28b) 
+ (a 2 +a 4 ) lU'MM 
+ ((ft -ft) |0»«1 K ,| + h.c) 

for the modes k and k? respectively. Notice that this 
formally corresponds to tracing "inside- out", that is, 
first (anti-)commuting operators towards the projector 
on the vacuum state before removing them, such that 

Tr m (bl\\0}(0 \\b m b n ) = | Mini ■ (29) 

We have now arrived at a point where we can make 
a general statement about the consistency conditions. 
Let us formulate this in the following theorem. 

Theorem 1. Given a density operator Qi t „ n G 
Ws{Fn) for n fermionic modes (labelled l,...,n) the 
consistency conditions (24) completely determine the 
reduced states on %s{Fm) for any m with 1 < m < n. 



Proof. This can be seen in the following way: for any 
matrix element 

lbl 1 ...bl\\0))(0\\b Vl ...b l/j (30) 

of an (n — l)-mode reduced state Qi c n -i) = 
T^nisi,. ..,«)) where 7 € C and the sets 

fx := {m,...,Hi}C {l,2,...,(n-l)} (31a) 

and v := {i^,...,^-} C {1,2,..., (n- 1)} (31b) 

label subsets of the mode operators for the (n — 1) 
modes, can have contributions from at most two matrix 
elements of Qi n > i- e - 

Tr n ( 7o 6t 1 ...6tj|0))((0|^ 1 ...6^), (32a) 

and Tr„ (71 b^ . . . b\ \\ )) (( || b n b Vx ...b Uj ). 

(32b) 

The composition of 7 into 71 £ C and 72 G C, i.e, 
7 = 7o±7i, is determined by the consistency conditions 
of Eq. (24). For every matrix element (30) with cor- 
responding partial trace contributions from (32) there 
exists a pair of hermitean operators 

O x {\,t)= b Xl ...b Xk b\ i ...b\ i (33a) 
+ b n ...b Tl b[ k ...b[ i , 

O p (X,t)= b Xl ...b Xk bl i ...bl l (33b) 
-ib n ...b Tl b[ k ■■■b{ i , 

with A := {Ai, . . . , A^} = p/v and t := {n, . . . ,n} = 
v/p,, that uniquely determine the relative sign of 71 
and 72. These operators are unique up to an overall 
multiplication with scalars. The tracing procedure can 
be repeated when any other of the (re — 1) remaining 
modes are traced over. Since the order of the partial 
traces is of no importance for the final reduced state, all 
reduced density operators are completely determined. 

□ 

Consequently, the reduced density matrices in the 
fermionic Fock space can be considered as proper den- 
sity operators, i.e., they are hermitean, normalized and 
their eigenvalues are well defined and non-negative. 
Moreover, since the eigenvalues are free of ambiguities, 
all functions of these eigenvalues, in particular all en- 
tropy measures for density operators, are well defined. 
Let us stress here that this analysis does not depend on 
any superselection rules that might be imposed in ad- 
dition. We will see how these enter the problem when 
mappings to qubits are attempted in Sec. VI. 

VI. ENTANGLEMENT OF FERMIONIC 
MODES 

We are now in a position to reconsider a measure 
of entanglement between fermionic modes. We can de- 
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fine the entanglement of formation E a p for fermionic 
systems with respect to a chosen bipartition A\B as 



E f(q) 



mm 

{Pn-|*„)) 



Pn £(\\ *„))), (34) 



in complete analogy to the usual definition [32]. Here 
the minimum is taken over all pure state ensembles 
W^n} that realize g according to Eq. (14) and £(|\&))) 
denotes the entropy of entanglement of the pure state 
I "J}}. Since the entropy of entanglement, e.g., using the 
von-Ncumann entropy, is a function of the eigenvalues 
of the reduced states Tr B ( || || ) or Tr A ( || || ) 

alone, we can conclude that this is a well defined quan- 
tity. As pointed out in Ref. [23], the minimization in 
Eq. (34) can be restricted to pure state decompositions 
that respect superselection rules. Since this restriction 
limits the set of states over which the minimization is 
carried out, the quantity without this restriction will be 
a lower bound to the "physical" entanglement of forma- 
tion. For two fermionic modes the minimization over 
all states that respect superselection rules can indeed 
be carried out, see Ref. [23]. However, in general this 
step will be problematic. 

Let us now turn to some operational entanglement mea- 
sure, in particular, let us investigate if and how the 
negativity N (see, e.g., Ref. [33]) and the concurrence 
C (see, e.g., Ref. [32]) can be computed to quantify 
fermionic mode entanglement. Both of these measures 
are operationally based on the tensor product structure 
of qubits. We will here define the negativity as 



(35) 



where are the eigenvalues of the partially trans- 
posed density matrix. The partial transposition is a 
map that is well defined only for basis vectors on a 
tensor product space. To employ this measure, let 
us therefore try to find a mapping 71", that is consis- 
tent with the conditions of Eq. (24). Starting with 
the two- mode state g KK i of Eq. (19) we are looking 
for a map tt that takes {|| )) , || 1 K )) , || 1 K > )) , || 1«)) || 1 K > ))} 
to {| 00), | 01), 1 10), 1 11)}, where \mn) = m ) ® 
| n ) G % K <8 H K >, such that 



(36a) 
(36b) 
(36c) 



A graphical representation of the consistency condition 
is shown in Fig. 1. It is quite simple to check that 
these requirements generally cannot be met, i.e., writ- 
ing g KK i of Eq. (19) as a matrix with respect to the 



Q >~ 




Qk i- 


-> T(e«) 


Qk' l- 


->■ 7r(e«') 



basis {|| )) , 1 1 K )) , I l K , )) , 1 1 K )) 1 1 K , ))} we get 



(ax Pi P2 Ps\ 

PI a 2 Pi Ps 

P*2 PI a 3 Pe 

\P* 3 PI PI aj 



(37) 



A mapping of the desired type should be obtained by 



Qkk 1 

g„ 



I Tr K 



FIG. 1. Graphical representation of consistent mapping: A 
mapping tt : g n> ir(g) from the space Hs{?2) to T-i K ®'H K i is 
considered to be consistent if it commutes with the partial 
trace operation. 

multiplying any number of rows and the corresponding 
columns by (—1) and considering the resulting matrix 
as the representation of ir{g) on H K ® H K >. The de- 
sired result should have a relative sign switch between 
Pi and Pq, while the signs in front of P2 and ft 5 should 
be the same. This clearly is not possible unless some of 
the coefficients vanish identically, e.g., by imposing su- 
perselection rules. For example, conservation of charge 
would require the coefficients /3i , /?2 5 /?s , fte, and, de- 
pending on the charge of the modes n and K 1 , either 
/?3 or P4 to vanish identically. In this way only inco- 
herent mixtures of pure states with different charge are 
allowed, but not coherent superpositions. 
We thus find that two fermionic modes can only be con- 
sistently represented as two qubits when charge supers- 
election is respected. In that case only one off-diagonal 
element can be non-zero and the sign of this element 
is insubstantial, i.e., it does not influence the reduced 
states or the value of any entanglement measure. In 
particular, the results for entanglement generation and 
degradation between two fermionic modes presented in 
Refs. [6, 8, 9] respect both charge superselection and 
the consistency conditions of Eq. (24). 
Let us return to the choice of entanglement measure 
for the permitted mappings to two qubits. We now re- 
strict the entanglement of formation E a p as defined in 
Eq. (34) to states that obey charge superselection, as 
suggested in Ref. [23]. As discussed earlier this means 
the usual entanglement of formation E Q p provides a 
lower bound to E Q p, i.e., 



< E, 



oF 



(38) 



For two qubits E p = E p(C) is a monotonically in- 
creasing function of the concurrence C. We propose an 
analogous functional dependence of E Q p = E p{C) on 
a parameter (7, that we call "fermionic concurrence" . 
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Evidently, the function C(g) is an entanglement mono- 
tone that is bounded from below by the usual concur- 
rence C. As shown in Ref. [34] the negativity N further 
provides a lower bound to the concurrence, i.e., in our 
convention of Eq. (35), 1M < C. Consequently, the 
negativity provides a lower bound to C, i.e., 

27V < C < C . (39) 

For two modes it is thus at least possible to compute 
lower bounds to entanglement measures explicitly. It 
was suggested in Ref. [20] that conventional entangle- 
ment measures overestimate the quantum correlations 
that can physically be extracted from fermionic sys- 
tems. The operations that can be performed on each 
single-mode subsystem are limited by (charge) supers- 
election as well. However, we conjecture that the inac- 
cessible entanglement between the fermionic modes can 
always be swapped to two (uncharged) bosonic modes 
for which the local bases can be chosen arbitrarily. 



VII. FERMIONIC ENTANGLEMENT 
BEYOND TWO MODES 

Finally, let us consider the entanglement between 
more than two fermionic modes. In principle, any mea- 
sure of entanglement that is based on entropies of the 
subsystems is well defined on the fermionic Fock space, 
as we have discussed. However, we would like to em- 
ploy operational measures. Let us therefore start by 
attempting a consistent mapping from three fermionic 
modes to three qubits, in analogy to the two-mode case 
in Sec. VI. For simplicity we assume that the modes k, 
k' and k" all have equal charge such that the most gen- 
eral mixed state of these modes can be written as 

Qn K 'n- - Mi II )) (( I + /ia I l K „ » {l K n I (40) 
+ II !«'»<( VI +M4 I V )) II V' K(l«" II « VII 

+ M5 |i K »((i K || + M6 IIUIIVMVIKUI 

+ M7 ||UIVK(VIK(VI 

+ /, 8 ||l K ))|l K ,))|V))((Vl((Vl((l K | 

+ (vi IV}) ((VII + »2 ||V ))((!« II 

+ ^3||V)}((l K | + ^|V}}|V}}((V|((l«l 
+ ^||M)|i K 4<(M<<VI 

+ ^ 6 IU|V>)<<M<(VI + h.c). 

The relevant consistency conditions to construct the 
three different reduced two-mode density matrices g KK > , 



g KK » and g K > K » are given by 

Tr((b{b n < + b\,b K ) eKK , K ») = 2Re(i*, + O , (41a) 

Tr((&tV + bl„b K )Q KK . K ,,) = 2Re(i/ 2 - i/ 5 ), (41b) 

Tr((blb K „ + bl,b K ,)g KK , K „) = 2Re(u 1 + i/ e ) . (41c) 

Again, the correct partial traces are obtained by tracing 
"inside-out", see Eq. (29). This is not a coincidence. 
The prescription for the partial trace to anti-commutc 
operators towards the projector of the vacuum state 
before eliminating them takes into account the number 
of anti-commutations occurring in computations of the 
expectation values of Eq. (24). A matrix representation 
of the three mode state Qkk'k" is given by 

M 0000000X 

fi 2 v\ v 2 

v\ (j, 3 i/ 3 

/14 U 4 I/ 5 . . 

0KK ' K " - v\ v% n 6 ' K > 

ul fia u 6 

z/| vl fj, 7 

Vooooooo^s/ 

Similar as before one can try to interpret Eq. (42) as 
a matrix representation of a three qubit state and ex- 
change the signs of the basis vectors in the three qubit 
state such that the consistency conditions of Eq. (41) 
are met, i.e., opposite signs in front of vi and while 
the signs in front of the pairs 1/3 , 1/4 and v\ , v§ are each 
the same. This is not possible, even though superse- 
lection rules are respected. This suggests that the su- 
perselection rules only coincidentally aid the fermionic 
qubit mapping for two modes. They simply force all 
the problematic coefficients to disappear. However, for 
more than two modes we find here that a mapping to a 
tensor product space cannot be performed consistently 
in general. Therefore, computing a measure like the 
negativity to determine the entanglement between more 
than two modes appears to be meaningless. Due to the 
lack of practical alternatives the minimization over all 
states consistent with charge superselection to find E f 
of Eq. (34) should be considered since the restriction of 
the set of permissable states could make this computa- 
tion feasible. 



VIII. CONCLUSION 

We have discussed the implementation of fermionic 
modes as fundamental objects for quantum informa- 
tion tasks. The foundation of this task is the rigorous 
construction of the notion of mode-subsystems in a 
fermionic Fock space. We have demonstrated, that 
this can be achieved despite the absence of simple 
tensor product structure. Our simple consistency con- 
ditions give a clear picture of this process, which can 
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operationally be easily executed by performing partial 
traces "inside-out". Thus we show that fermionic 
mode entanglement, quantified by the (fermionic) 
entanglement of formation, is indeed a well defined 
concept, free of any ambiguities. 

However, problems arise when mappings from the 
fermionic Fock space to qubit spaces arc attempted. 
We have explicitly demonstrated in two examples, 
for two and three modes, that such mappings cannot 
generally succeed. Only in the limited case where only 
two modes are considered and the quantum states 
obey charge superselection can one meaningfully speak 
of an equivalence between the two fermionic modes 
and two qubits. In this case the application of tools 
such as the negativity or concurrence is justified. We 
have argued that these measures will at least provide 
a lower bound to genuine measures of fermionic mode 
entanglement. 

Nonetheless, open questions remain. In particular, it is 
not clear if any operational measures exist for situation 
beyond two qubits. In Ref. [10] witnesses for genuine 
multipartite entanglement are employed, which are 



completely compatible with the framework we have 
presented here, but these witnesses can only provide 
lower bounds on entropic entanglement measures. 
Finally, we have conjectured that the entanglement 
in fermionic modes is accessible even in spite of su- 
perselection rules that restrict the possible operations 
performed on single modes by means of entanglement 
swapping. The investigation of this question, while 
beyond the scope of this article, will certainly be of 
future interest. 
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